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The ability to assemble nano- and micro- sized colloidal components into highly ordered configurations is 
often cited as the basis for developing advanced materials. However, the dynamics of stochastic grain 
boundary formation and motion have not been quantified, which limits the ability to control and anneal 
polycrystallinity in colloidal based materials. Here we use optical microscopy, Brownian Dynamic 
simulations, and a new dynamic analysis to study grain boundary motion in quasi-2D colloidal bicrystals 
formed within inhomogeneous AC electric fields. We introduce "low-dimensional" models using reaction 
coordinates for condensation and global order that capture first passage times between critical 
configurations at each applied voltage. The resulting models reveal that equal sized domains at a maximum 
misorientation angle show relaxation dominated by friction limited grain boundary diffusion; and in 
contrast, asymmetrically sized domains with less misorientation display much faster grain boundary 
migration due to significant thermodynamic driving forces. By quantifying such dynamics vs. compression 
(voltage), kinetic bottlenecks associated with slow grain boundary relaxation are understood, which can be 
used to guide the temporal assembly of defect-free single domain colloidal crystals. 

Photonic and meta- materials provide examples where ordered particles on length scales comparable to 
electromagnetic wavelengths produce exotic emergent properties \ Colloidal crystallization provides a 
potential route to self- assemble such materials via processes amenable to scalable manufacturing; however, 
robust schemes have yet to be identified to obtain the necessary low defect densities. From a broader perspective, 
obtaining perfect crystals on any length scale remains more art than science {e.g., atoms, molecules, macromo- 
lecules)^'\ To design robust crystal growth, recrystallization, and annealing schemes to minimize defects, it is 
necessary to know the basic mechanisms of defect formation and motion. 

Although three-dimensional configurations of complex particles are the ultimate goal of self-assembly 
schemes^'^, here we investigate a relatively simple problem that is still not well understood: how grain boundaries 
form and move during quasi two-dimensional (2D) crystallization of spherical colloids. Quasi-2D colloidal 
crystals have been used in studies of melting^, nucleation^, point defect diffusion^, and grain boundary fluctua- 
tions^. Other relevant studies include 2D analyses of colloidal crystals to investigate impurity mediated growth^°, 
particle motion within grain boundaries^ \ and grain boundary pre-melting^^. Despite these extensive studies, 
fundamental understanding of grain boundary formation and motion remains rudimentary^^ in all but the most 
model simulation studies^^. 2D crystals are also of interest based on their relevance to thin fllms^^ bubble rafts^^, 
and graphene^^. Understanding how grain boundary motion enables relaxation of multi-domain crystals into 
defect- free crystals is therefore scientifically and technologically interesting. 

Results 

In this work, we employ real-time microscopy to observe grain boundary formation and motion in a quasi-2D 
colloidal crystal containing —200 colloidal particles (Fig. la. Supplementary Videos 1-4). Aqueous ~3 micron 
Si02 charged colloids crystallize in a quadrupole electrode in MHz AC electric fields (see Methods, 
Supplementary Information) ^^ where field-mediated compression of induced dipoles is balanced by the quasi- 
2D colloid osmotic pressure (i.e., effective hard disk). The electric field amplitude (i.e., applied voltage) is 
effectively a surrogate for pressure and hence acts as a global thermodynamic variable that determines the relative 
free energy of all particle configurations for fixed voltage, number, and temperature. In the following, we report 
non-dimensional voltages, where V* = 1 is the voltage required for N particles to produce a hexagonal close 
packed crystal with hexagonal morphology (as demonstrated in previous work with agreement between micro- 
scopy experiments, MC simulations, and perturbation theory^^, see Supplemental Information for additional 
details). 
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Figure 1 | Reaction coordinates for grain boundary formation and motion can be computed from image analysis of optical microscopy images of 
electric field mediated colloidal crystallization, (a) Raw images show representative configurations of 210 ~ 3 \im silica colloids in fluid, bicrystal, and 
single crystal configurations observed over —7.5 min following a quench (step voltage change) to V* = 0.57 (see main text and Supplementary Methods 
for definition). Snapshot times are shown by solid black lines in Fig. le. Computed reaction coordinates shown by colored particle centers on 8 -bit 
intensity scale for the (b) radius of gyration, i^g, (c) local hexagonal order, (Ce), (d) global hexagonal order, i//6, and time dependent traces for (e) a single 
voltage quench and (f) ten consecutive cycles (numbered vertical black lines indicate representative images included in Supplementary Fig. 2). 

Our previous characterization of electric field mediated colloidal icles^^. Based on these findings, we set AT = 210 in Fig. 1, which 
interactions and assembly has yielded /cT-scale potentials^°'^\ feed- routinely forms bicrystals (i.e., 1 grain boundary between 2 domains) 
back control over system size^^, and conditions to crystallize AT part- in contrast to single domains in smaller systems and >2 domains in 
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larger systems. As shown in Fig. 1, step-quenches to = 0.57 cause 
an initially dilute fluid phase to first rapidly condense, then form 
grain boundaries via coalescence of local domains, and finally display 
grain boundary motion as bicrystals relax to single crystals. 
Reversibility allows repeated quenches between fluid and crystal 
states to probe the stochastic dynamics of grain boundary formation, 
diffusion {i.e., random motion), and migration {i.e., drift). We first 
focus on the V* = 0.57 case, but later investigate step-changes to 
both lower and higher values of V^. 

To interpret and model these measurements, we aim to develop a 
"low- dimensional model" that quantitatively captures the observed 
dynamics using "reaction coordinates"^^ (rather than enumerating 
all 2Artranslational degrees of freedom) The use of the term "low- 
dimensional" here does not refer to the Euclidian spatial dimension 
{i.e., X and y particle center coordinates), which is effectively quasi- 
two dimensional, but rather, "dimensionality" refers to the number 
of reaction coordinates necessary to capture the stochastic dynamics 
of grain boundary formation. It is "low dimensional" because we 
expect the number of reaction coordinates to be considerably less 
than the 2 X (AT = 210) = 420 dimensions that would be necessary to 
uniquely specify all possible two-dimensional configurations of 210 
particles. Although the term "low- dimensional" may be unfamiliar 
to some readers, other synonymous terms such as "coarse-grained 
models" could also be confusing based on an unfortunate overlap of 
terms relevant to the application of interest in this work. 

Candidate reaction coordinates are computed from particle cen- 
ters and used to color code images (Figs. Ib-f) including: the radius 
of gyration, R^^, to capture condensation from fluid to crystal states, 
average local hexagonal order, {C^^^, to capture the onset of crystal- 
lization, and global hexagonal order, xjje^, to capture the degree of 
polycrystallinity. and {C^ are normalized by their AT particle single 
crystal values (see Supplementary Methods), so that {C^ goes from 
0-1 for fluids to complete locally ordered states, and R^ decreases 
from arbitrarily high numbers to 1 for complete condensation^^. is 
0 for fluids and 1 for single-domain crystals like {C^, but in contrast, 
depends strongly on relative domain size and misorientation {e.g., i/^6 
= 0 for bicrystals of identically sized 111 domains with 30° 
misorientation). 

The reaction coordinates trajectories following a single voltage 
quench (Fig. le) show initially decreasing R^ and increasing {C^ 
occur in unison {i.e., Rf^ ~ (Q))^ which demonstrates a close coup- 
ling between condensation and local ordering. As R^ and {C^ plat- 
eau, indicating an overall condensed conflguration with all particles 
contained in locally crystalline domains, i/^6 ~ 0» indicating a bicrys- 
tal. For t > 100 s, i^g and {C^ remain essentially unchanged while i/^6 
rises from 0 to —0.8, which clearly corresponds to grain boundary 
motion from the bicrystal interior to the periphery where it vanishes 
(Fig. Id, Supplementary Video 1). 

Ten successive voltage quenches from initial fluid states (Fig. If) 
demonstrate the stochastic nature of the grain boundary dynamics, 
which is expected from the underlying probabilistic colloidal motion. 
While i^g and {C^ reveal condensation and local order emerge in a 
similar manner for each cycle, the xj/^ trajectories can be categorized 
into several cases: (1) xj/^ tracks {C^ indicating the simultaneous 
emergence of local and global order, (2) becomes localized at 
intermediate values for varying time periods before again increasing, 

(3) i/^6 becomes arrested for the duration of the observation time, and 

(4) in one case, i/^6 initially increases but then vanishes. The stochastic 
nature of the observed grain boundary dynamics apparent in an 
ensemble of trajectories is an important aspect to capture in a quant- 
itative model. 

To develop a low- dimensional model of grain boundary formation 
and motion, it is necessary to determine the number and type of 
reaction coordinates. For example, it could be speculated that i/^6 is 
all that is required to track grain boundaries since it visually tracks 
polycrystallinity in Fig. 1. However, simply tracking does not 



capture how parallel processes of local condensation and crystalliza- 
tion determine the formation and motion of grains of different sizes, 
shapes, and orientations, which ultimately determine the mechan- 
isms of polycrystals relaxing to single crystals. In addition, 1/^5 does 
not uniquely identify some conflgurations; for example, fluid con- 
flgurations and maximally misaligned bicrystals both have {1/^ = 0. 

To illustrate how reaction coordinate pairs capture additional 
information, experimental conflgurations are colored using two 
coordinates (Figs. 2a, b) and 2D trajectories of {ij/^, (C^)) (Fig. 2c) 
or {ij/^, i^g) (Fig. 2d) vs. time, where time is indicated by a 256-color 
scale. We do not plot {R^, (C^)) since these coordinates are highly 
correlated in Fig. 1, although such a pair could be useful to distin- 
guish condensed amorphous microstructures {i.e., glasses/gels). Such 
2D plots immediately address one issue; plotting ij/^ against either R^ 
or (Ce) distinguishes fluid and bicrystal conflgurations {e.g., high Rg, 
low 1/^6 vs. low Rg, low ij/^). These trajectories also show how local 
ordering {i.e., increasing (C^)) and condensation {i.e., decreasing R^) 
influence the emergence of polycrystallinity and subsequent grain 
boundary motion. Trajectories starting at lower ij/^ tend to become 
localized at lower ij/^ after condensation, whereas trajectories with 
initially higher ij/^ tend to rapidly form single crystals. 

Although tools exist to identify the minimum dimensionality {i.e., 
number of reaction coordinates) from observed dynamics {e.g., dif- 
fusion mapping^^"^^), such methods are not currently able to predict 
physically meaningful reaction coordinates^^. Ultimately, the num- 
ber and types of reaction coordinates can be determined empirically 
by finding what is necessary to produce a quantitative stochastic 
dynamic model. Because the measured trajectories display both drift 
and diffusion along reaction coordinates, which appear to be 
mediated by free energy gradients {i.e., driving forces) and fluctua- 
tions/friction {i.e., randomness/resistance), it is assumed that such 
processes can be captured by a low- dimensional Smoluchowski 
equation given by^°. 



dp{x,t) 
dt 



^V-D(x)- 



V+^VW(x) 
Kb! 



p{x,t) 



(1) 



where p{x, t) is the probability density of finding the system at coord- 
inate x at time, t, W{x) is the free energy landscape, kT is thermal 
energy, and D(x) is the diffusivity landscape, which can be related in 
the usual way to mobility, m(x), and friction, ^(x), landscapes as D(x) 
= m(x)A;T = l^{x)~^kT. The variable x is a vector of reaction coordi- 
nates where the number of coordinates is the model "dimension." 
The Smoluchowski equation is a special case of the Fokker- Planck 
equation that also satisfies the fluctuation-dissipation theorem and 
leads to Boltzmann sampling at equilibrium {i.e.,p{x) = exp[- W(x)/ 
^23,30 more descriptive terms, W{x) is the free energy change 
associated with moving from one particle configuration to another 
(or one reaction coordinate to another; {e.g., Xi = (i/^6,i» ^g,i) to ^2 = 
{^6,2^ ^g,2))> and D(x) captures the associated configuration depend- 
ent changes in diffusion and friction {i.e., fluctuations and 
dissipation). 

To obtain W{x) and D(x) in Eq. (1) from trajectories like those in 
Figs. 1 and 2, we analyzed Brownian Dynamic (BD) simulations that 
were matched to experiments (by capturing all equilibrium and 
dynamic properties of the quadrupole experiment on the particle 
scale see Supplementary Methods). This approach was used 
because statistics on the particle scale are easily obtained to match 
experiments and simulations, but BD simulations are better suited to 
generating large numbers of grain boundary trajectories {e.g., each 
experimental grain boundary trajectory in Fig. 1 is acquired for 
~10 min.). 

To provide more details of the matching process, inverse Monte 
Carlo was used to obtain interaction potentials that capture all equi- 
librium properties {i.e., radial distribution functions), which were 
then employed in BD simulations that captured all dynamic prop- 
erties {i.e., particle scale diffusion, reaction coordinate trajectories). It 
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Figure 2 | "Two dimensional" trajectories (i.e., two reaction coordinates) capture coalescence of local domains during fast condensation processes to 
produce bicrystals that relax to single crystals via grain boundary motion over a broad range of timescales. Microscopy images from Figs, la with 
particle centers colored using RGB (Red/Green/Blue) mixing rules for colors represented by two reaction coordinates including (a) ips and Cg and (b) ij/s 
and i^g to visualize how global order emerges from local order and during condensation. Ten trajectories following quenches to V* = 0.57 with time 
represented by a 256-color scale (inset scale bar) for (c) (i//6> (Q)) and (d) (i//6> reaction coordinate pairs with inset plots of single trajectories 
from Fig. le. 



is important to note that it was necessary to include concentration 
dependent, but not configuration dependent, hydrodynamic inter- 
actions in the BD simulations to produce agreement with experi- 
ments. In particular, Stokesian Dynamic simulations^^ were used to 
show that approximate local, concentration dependent hydrodyn- 
amic interactions were sufficient to match experiments without the 
need to consider multi-body, configuration dependent hydrodyn- 
amic interactions. This finding indicates local microstructure or 
morphology do not obviously influence the hydrodynamic interac- 
tions involved in the colloidal crystal grain boundary dynamics 
observed in our work. 

Statistical methods reported in the literature^^ and further 
developed by us for application to colloidal assembly^^, were used 
to analyze large numbers of BD simulated trajectories to construct 
W{x) and D(x) (see more details in Supplementary Methods and our 
previous work^^'^^). In brief, the displacement and mean squared 
displacement of reaction coordinate vs. time trajectories can be used 
to measure drift and diffusion at each value of x, which ultimately 
yield W(x) and D(x).To assess the quantitative accuracy of candidate 
low-dimensional dynamic models, we compared first passage time 
distributions for ensembles of trajectories between different starting 
and ending states from particle-scale BD simulations and low- 



dimensional Langevin dynamic (LDLD) simulations. The LDLD 
simulations are based on a Langevin equation given as, 

x(t + At) = x{t)-D[x{t)]'{kTy\w[x{t)]At 

(2) 

+ V-D[x(0] + [2D[x(0]At]^/'r(0 

where the coefficients are the same as in Eq. (1), At is the integration 
time step, and T(t) is a noise variable. A successful LDLD model will 
accurately reproduce first passage time distributions obtained from 
high 2 AT- dimensional particle- scale BD simulations, and hence the 
experimental trajectories that were quantitatively matched to the BD 
simulations. 

The only quantitatively accurate model to emerge from candidate 
ID and 2D models was one based on (i/^g, R^) (so higher dimensional 
models were not considered). Fig. 3a shows the W{x) that quantifies 
the relative free energy of every configuration and the free energy 
gradients that drive motion. The global free energy minimum at x = 
(i/^6» ^g) ~ (0-8, 1.14) indicates a single domain crystal (with a thin 
fluid envelope at its periphery) is the thermodynamically favored 
configuration (image VI in Fig. 3c). 

The diffusivity landscape, D(x), in Eqs. (1) and (2) is a 2 X 2 
tensor; it has diagonal components that capture how friction/mobil- 
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Figure 3 | Free energy and friction landscapes obtained by fitting trajectories to Smoluchowski equation (equation 1) that quantitatively capture grain 
boundary diffusion and migration vs. relative domain sizes and misorientation angles, (a) W(i//6, Rg)/kT with inset scale bar and two trajectories 
obtained from the experiments in Figs. 1 and 2. (b) T>^^^J{kT-s) with inset scale bar with the same trajectories as in (a), (c) Representative configurations 
from microscopy images and simulated renderings for coordinates marked along trajectories in (a) and (b). Labeled misorientation angles shown by red 
lines and particle centers are colored according to the RGB composite convention in Fig. 2. First passage time distributions for BD (red) and LDLD (black) 
trajectories projected onto (d) the i/^e-axis between 0.47-0.66 (#,•), 0.38-0.56 (A,A), 0.28-0.47 (■,■), 0.19-0.38 (♦,♦) and (e) i^g-axis between 
L27-1.25 (#,#)a.25-1.23(A,A)a.23-1.21 (■,B)a.21-l. 19 (♦,^)a.21-l. 19 (#,#)a. 19-1. 17 (TT), and (f) for trajectories between a sink 
= 0.8 i^g = 1.14) {i.e., global minimum) and sources at (1/^5 = 0.38 = 1.15) (#,•) and (1/^6 = 0.65 = 1.16) (A,A). 



ity in each reaction coordinate mediate drift due to free energy gra- 
dients in the same coordinate (D^6\|/6> ^RgRg), and cross-terms that 
capture how friction/mobility mediate drift due to free energy gra- 
dients in orthogonal coordinates (D^6Rg» ^Kg^e)- Fig- 3b shows 



Dx|/6\|/6/^^j which captures how friction mediates drift and diffusion 
in 1/^6 due to free energy gradients in The other components of 
D(x) are reported in Supplementary Fig. 5. DRgRg shares similar 
features with D^6\|/6- The cross-terms (D^6Rg» DRg^e) indicate a weak 
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coupling between driving forces and drift/diffusion (i.e., friction 
increases as free energy increases) along orthogonal coordinates 
for il/^>0.5, although this is relatively minor compared to the diag- 
onal terms. 

Representative configurations (Fig. 3c) and first passage time dis- 
tributions (Fig. 3d-f) show the resulting low dimensional il/^, 
model quantitatively captures the measured grain boundary 
dynamics. The agreement between the BD and LDLD simulations 
is excellent (Fig. 3d-f), demonstrating that Eqs. (1) and (2) with the 
W(i/^6> ^g) and D(i/^6> ^g) in Fig. 3 provide accurate low dimensional 
dynamic models of the experiments in Figs. 1 and 2. 

Discussion 

To aid discussion of how features on W(i/^6> ^g) and D(i/^6» ^g) are 
connected to microscopic mechanisms, two limiting trajectories are 
shown on these landscapes; one where a grain boundary forms 
between two domains and does not move out of the crystal on the 
— 10 min observation time [trajectory 1 (Tl): I-II-IV-V], and one 
where two domains form but the grain boundary quickly moves to 
the crystal periphery in < 1 min to produce a single domain crystal 
[trajectory (T2): I-III-VI]. Both trajectories are consistent with 
expectations by showing drift (i.e., migration) along free energy gra- 
dients and superimposed stochastic motion (i.e., diffusion) that is 
most evident where free energy gradients are minimal. The dramatic 
difference between these two trajectories is most evident as they 
approach the global minimum (Fig. 3f) (as quantified between 
"sources" at III and IV and a "sink" at VI). An order of magnitude 
difference is observed in the most probable first passage times with 
-20 s for III-VI and -200 s for IV- VI; the latter distribution also 
shows a much longer asymmetric tail with some trajectories taking 
> 1,000 s to traverse the W(i/^6» ^g) plateau at low R^. 

The Tl trajectory corresponds to rapid condensation along a steep 
free energy gradient where two locally ordered domains coalesce into 
a bicrystal with a near maximum 30° misorientation angle. At point 
IV on Tl, the grain boundary randomly diffuses on a free energy 
plateau with a minimal free-energy gradient (i.e., driving force) to 
drive migration of the grain boundary to the crystal edge; the tra- 
jectory is locaUzed between IV-V for —10 min. In addition to van- 
ishing free energy gradients, the friction in the vicinity of IV-V is 
increased — 6X compared to uncondensed states. Multi-body hydro - 
dynamic interactions (i.e., near-field lubrication and far-field flow 
within the particle structure^^) increase particle-scale friction during 
condensation, which is consistent with the increased friction for 
trajectories at low Rg on D(\l/^, Rg) (Fig. 3b). Although first passage 
times vary linearly with frictional changes compared to an exponen- 
tial dependence for free energy changes^^'^^ diffusion mediated by 
friction is the rate determining process in the presence of vanishing 
free energy gradients (plateau of Fig. 3a). 

In contrast to Tl, the T2 trajectory is initiated with higher global 
order before moving down the free energy gradient and rapidly 
continuing towards the global free energy minimum single crystal. 
Although T2 passes close to Tl, low friction at high R^ allows suf- 
ficient diffusion towards higher ij/^. As a result, T2 bypasses the free 
energy plateau region at low Rg to avoid slow diffusion like Tl, and 
instead shows much faster grain boundary migration. Friction uni- 
formly increases with decreasing Rg (due to hydrodynamic interac- 
tions), and has almost no dependence on ij/^, so there is no path of 
least resistance on D(\l/^, Rg). In short, the fastest trajectories are ones 
that bypass the free energy plateau region. 

The microscopic mechanisms associated with these different tra- 
jectories can be understood from the images/renderings (Fig 3c) and 
the physical meaning of the reaction coordinates. Rg clearly captures 
condensation as shown by large free-energy gradients on W(i/^6» ^g) 
and increasing resistance to configurational changes on D(\l/^, Rg) as 
the result of multi-body hydrodynamic interactions. At this point, we 
can speculate why (C^) was not part of a successful dynamic model; it 



is an indirect measure of condensation, and thus not as good as Rg, 
and it is the emergence of global order, captured by ij/^, that is most 
important to track grain boundaries. 

The i/^6 dependence of W(\l/6, Rg) indicates that domains coales- 
cing with minimal misorientation produce higher global order from 
the outset, which also translates to faster grain boundary migration 
(via free-energy gradients) from the crystal interior to the periphery. 
Practically, low misorientation angles produce smaller energy bar- 
riers to particle-scale motion (in full 2 AT- dimensional particle-scale 
space) within (e.g., string-like motion) and across (e.g., cooperative 
motion) grain boundaries'^. In contrast, domains that coalesce near 
the maximum 30° misorientation display low initial global order that 
translates into slow grain boundary diffusion on a W(i/^6» ^g) plateau. 
Such bicrystals represent an unstable equilibrium where the energy 
(e.g., energy/atom, interfacial energies) of the two sides balance, 
however, the lower free energy state single crystal (—10 kT) is 
achieved by fluctuations that eventually allow one grain to increase 
at the expense of the other grain decreasing. 

Because i/^6 does not resolve different combinations of domain size 
and misorientation, but accurately captures the dynamics, it appears 
that all such configurations relax in an indistinguishable manner. In 
particular, greater misorientations between dissimilar sized grains 
produces relaxation rates equivalent to cases where domains have 
less misorientation but are of similar size. This finding shows how 
using i/^6 as a reaction coordinate indicates an aspect of grain bound- 
ary motion that would not be easily discovered from tracking par- 
ticle-scale motion alone. Because grain boundary motion involves 
many particles rearranging in a cooperative fashion based on relative 
domain sizes and orientations, it is useful to have a global parameter 
that captures configurational changes of the entire particle ensemble, 
and therefore naturally captures cooperative phenomena. Although 
i/^6 should be a good reaction coordinate to monitor relaxation of 
bicrystals of any size (i.e., distinguish two-domain from single- 
domain crystals), it may be unsuitable for many-domain crystals 
because it would remain near zero and miss most of the structure 
evolution until only a few domains remained. In many- domain crys- 
tals, it might be more appropriate to capture the evolution of poly- 
crystallinity using a different reaction coordinate, like average 
domain size or average misorientation angle across grain boundaries. 

The approaches used to measure and model grain boundary 
motion for V* = 0.57 in Figs. 1-3 can be applied at other V*. 
Using the same BD simulations and non- equilibrium analyses, 
W(il/6, Rg) were constructed in the range = 0.31-0.69 (Fig. 4a- 
d) along with representative trajectories and global minimum con- 
figurations. At the lowest V^, particles are weakly confined in con- 
centrated fluid configurations without crystal grains or boundaries, 
which produces a relatively featureless W(il/6, Rg). At the highest V^, 
W(il/6, Rg) is qualitatively similar to Fig. 3a, but the plateau stretches 
from very low to high ij/^ with an even shallower gradient and deeper 
global minimum at i/^6 ^ 1- At intermediate V*, the W(\l/6, Rg) show 
a continuously shifting global minimum towards lower Rg and higher 
ij/e and a stretching plateau corresponding to slower grain boundary 
migration. D(x) vs. are not reported here, but the general trend 
with increasing voltage is a decreasing magnitude (i.e., decreasing 
mobility, increasing friction), which is consistent with more con- 
densed configurations hindering particle motion and hence motion 
along both reaction coordinates. In short, with increasing compres- 
sion, once grain boundaries form, they experience slower migration 
and diffusion, and the single perfect crystal clearly emerges as the 
global free energy minimum configuration. 

The results show voltage cannot simply be increased to increase 
order, since in general this produces increasingly arrested polycrys- 
talline states. It is interesting to consider how grain boundary forma- 
tion and motion might be manipulated by "switching" between 
W(il/^, Rg) at different V*. If is increased very slowly, it would 
be possible to remain in the global free energy minimum configuration 
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Figure 4 | Free energy landscapes based on i//^,^ reaction coordinate 
pair capture dynamics at all V* to provide quantitative models of grain 
boundary formation and motion. W(iA6> R^)lkTdX (a) V* = 0.31, (b) V* = 
0.44, (c) V* = 0.57, and (d) V* = 0.69 with inset renderings of global 
minimum configuration and representative trajectories from BD 
simulations that were matched to the experiments (see Supplementary 
Methods). The free energy scale is indicated by the inset in (a). 

at all voltages; this is the thermodynamic equilibrium limit and a 
known strategy to make single crystals. However, faster schemes are 
generally desirable. The optimal control policy^^ to achieve a single 
perfect crystal in minimal time could employ the quantitative non- 
equilibrium dynamic models in Eqs. (1) and (2) to switch between 
W(i/^6» ^g) at different in an automated, informed manner using 
feedback controP^ In particular, monitoring reaction coordinates in 
real-time could identify slowly relaxing polycrystalline configura- 
tions {i.e., due to vanishing free-energy gradients and high friction), 
and then V* could be tuned to "land" on another W(i/^6» ^g) at the 
same coordinates where faster relaxation occurs. 

In summary, we report agreement between optical microscopy 
measurements, Brownian Dynamic simulations, and low-dimensional 



models of stochastic grain boundary formation and motion in quasi- 
2D colloidal bicrystals. Our results show that two reaction coordi- 
nates, one for condensation and one for global order, are sufficient to 
quantitatively capture first passage times between critical configura- 
tions at each applied voltage. Free energy and diffusivity landscapes 
show that the relative misorientation angles and domain sizes 
formed during condensation determine the subsequent grain bound- 
ary motion. Bicrystals with similar sized domains and a near 30° 
maximum misorientation angle relax via slow grain boundary dif- 
fusion mediated by high friction and vanishing free energy gradients, 
whereas bicrystals with asymmetrically sized and/or less misoriented 
domains relax via much faster grain boundary migration due to 
greater thermodynamic driving forces. By quantifying such 
dynamics as a function of voltage, ongoing work is developing 
optimal control algorithms to dynamically tune voltages to avoid 
kinetic bottlenecks associated with slow grain boundary dynamics. 

Future work will extend the modeling approaches developed in 
this work to understand other mechanisms and defect types involved 
in colloidal crystal formation/relaxation {e.g., nucleation, growth, 
point defects, multiple domains, etc.), which will require identifica- 
tion of the minimal number and type of order parameters to accur- 
ately capture relevant dynamic processes. Additionally, the role of 
multi-body, configuration dependent hydrodynamic interactions on 
diffusivity (friction) landscapes will be explored in attractive particle 
systems, where changes in near-field interactions during clusterings^ 
and percolation^^ have previously been shown to have non-trivial 
dynamic signatures. 

Methods 

Coplanar gold thin film quadmpole electrodes were patterned on glass microscope 
coverslips by spin coating photoresist and physical vapor deposition of a 15 nm 
chromium layer and a 35 nm gold layer. Nominal 3.13 |im diameter Si02 colloids 
with —50 mV zeta potentials were fractionated in Dl water and centrifuged/redis- 
persed five times in 0.1 mM NaOH. PDMS o-rings were coated with vacuum grease 
and sealed between a coverslip with the patterned quadrupole electrode before it was 
connected in series with a function generator. Microscopy was performed on an 
inverted optical microscope with a 63X objective and a 12-bit CCD camera that 
captured 336 pixel X 256 pixel (81 |im X 62 |im) digital images at rate of 8 frames/s. 
Video capture and image manipulation were performed using algorithms in 
MATLAB. 

BD simulations in the canonical ensemble were performed for 210 colloidal part- 
icles at constant voltage using numerical methods described in previous 
papers^^'^^'^°"'*^. A 0.1 ms time step was used for at least 2 X 10^ steps, and reaction 
coordinates were stored every 1250 steps for subsequent analysis. Particles in simu- 
lations were confined within 2D planes. Inverse Monte Carlo methods (including 
image resolution limiting effects)^^""^^ were used to match measured and simulated 
radial distribution functions to determine parameters in interactions potentials. The 
diffusivity was matched by comparing measured and simulated mean square dis- 
placements. Parameters used in the BD simulations are reported in the 
Supplementary Information. 
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